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' — ' ABSTRACT 
> 

0^ We present deep emission-hne imaging taken with the new SOAR Optical 

^ Imaging Camera of the brightest cluster galaxy (BCG) in the nearby {z = 0.035) 

^ — J X-ray cluster of galaxies 2A0335+096. We also present our analysis of addi- 

^ tional, multi-wavelength observations for the BCG, including long-slit optical 

f — spectroscopy, archival VLA radio data, Chandra X-ray imaging, and XMM UV- 

^ imaging. Cluster 2A0335-I-096 is a bright, cool-core X-ray cluster, once known 

as a cooling flow. Within the highly disturbed core revealed by Chandra X- 
^ ray observations, 2A0335-I-096 hosts a luminous and highly structured optical 

a emission-line system, spanning the brightest cluster galaxy (BCG) and its com- 

panion. Wc confirm that the redshift of the companion is within 100 km s~^ of 
the BCG and has certainly interacted with the BCG, and is likely bound to it. 
The comparison of optical and radio images shows curved filaments in Ha emis- 
sion surrounding the newly resolved radio source. The velocity structure of the 
emission-line bar between the BCG nucleus and the companion galaxy provides 
strong evidence for an interaction between the BCG and its northeast companion 
in the last ~ 50 miUion years. The age of the radio source is similar to the inter- 
action time, so this interaction may have provoked an episode of radio activity. 
We estimate a star formation rate of > 7 M© yr~^ from the Ha and archival UV 
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data. This rate is similar to, but somewhat lower than, the revised X-ray cooling 
rate of 10 — 30 Mq yr~^ in the vicinity of the BCG, estimated from XMM spectra 
by Peterson et al. (2003). The Ha nebula is limited to a region of high X-ray 
surface brightness and cool X-ray temperatures. However, the detailed structures 
of Ha and X-ray gas differ. The peak of the X-ray surface brightness is not the 
peak of Ha emission, nor does it lie in the BCG. The estimated age of the ra- 
dio lobes and their interaction with the optical emission-line gas, the estimated 
timescale for depletion and accumulation of cold gas, and the dynamical time in 
the system are all similar, suggesting a common trigger mechanism. 

Subject headings: galaxies:clusters:general — galaxies:clusters:individual (2A0335+096) 
— cooling flows 



Introduction 



In the last few years, remarkable progress in our understanding of the properties of hot 
gas in the centers of clusters of galaxies has been made possible by the Chandra X-ray and 
XMM (X-ray Multi-Mirror) observatories. Stringent constraints on emission lines from the 
appropriate ionization species of iron and oxygen in high-resolution X-ray spectra from XMM 

but indicate the 



show little evidence for simple, high M cooling flows (Peterson et al. 



2003 



presence of gas ranging over a factor of two in X-ray temperature. High resolution imaging 
(to ~ 1") by Chandra has revealed cavities in the intracluster medium (ICM), inflated by 
radio plasma ( [McNamara et al.||2001| |2000| [Birzan et aL]|2004D . The lack of hot rims around 
the bubbles suggested that the inflation has been gentle. Chandra spectroscopic imaging has 
also revealed X-ray fronts, which are a likely to be the manifestation of cooler X-ray blobs 



passing through hotter gas (e.g. Markevitch et al. 2000). Very deep Chandra observations 



have also revealed shocks emanating from the center of at least two clusters (McNamara et al. 



2005 Nulsen et al. 2005). Within 10 kpc of the core, structures in the hot X-ray emitting 
plasma often trace structures seen in optical emission-line images (e.g., Fabian et al. 2001 



McNamara et al. 2001; Sparks et al. 2004). 



AGN feedback may be the key to stabilizing cool cores in X-ray clusters. Such expla- 
nations have become increasingly popular to account for the high frequency of clusters that 
have a large amount of gas with a short radiative cooling time, multi-temperature gas, and 
non-radial spatial structures in the cluster core. The earliest hydrodynamic feedback models 



were produced to explain the heating and cooling of X-ray gas in elliptical galaxies (Binney 



& Tabor 1995 Tabor & Binney 1993). With the Chandra discoveries of bubbles in many 



nearby clusters, these models for elliptical galaxies evolved into feedback models for the ICM 
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massive black holes in the universe, at the centers of clusters 
uce winds, jets, and explosions that alter the ICM entropy and 


(Lauer et al.||200( 


3), may proc 



prevent catastrophic cooling from occurring in clusters. 

AGN feedback also appears to be an increasingly popular and important ingredient in 
galaxy formation scenarios. Mergers may couple the growth of galaxies and their central 
black holes, thus providing a natural explanation for the unexpectedly tight relationship. 



discovered by Gebhardt et al. (2000) and Ferrarese & Merritt (2000) between the masses of 
those black holes and their host's bulge luminosity and the velocity dispersion (e.g. Haehnelt 



& Kauffmann 2000). Mergers alone can not reproduce the exponential cutoff of the galaxy 



luminosity function, but mergers with AGN feedback might shut off star formation in very 



large galaxies (e.g. Best et al. 2006 Croton et al. 2006). The AGN heating we now ob- 
serve in low redshift systems may therefore give us insight into how AGN heating (and 
feedback) might have worked in the distant past. Semi-analytic studies currently make gen- 
eral assumptions about the efficiency of AGN feedback, tuned to suppress the formation 
of over-luminous, blue galaxies at the centers of clusters and to reproduce the exponential 
cut-off of the galaxy luminosity function at the high end (e.g. Croton et al. 2006 Best et al. 



2006). But these studies, because of the need for computational efficiency, must omit the 



details of how the feedback works. State-of-the-art hydrodynamic models of individual radio 



sources generally neglect the complexity of the ambient ICM (e.g. Reynolds et al. 2005 De 



Young & Jones 2006). It is necessary, therefore, to test the assumptions going into these 
important models with detailed observations of individual systems, in order to further the 
goal of making more detailed models of such systems, including triggering and termination 
of AGN heating. 

We describe here multi-wavelength observations of the X-ray cluster 2A0335-I-096 that 
illustrate how complex these interactions can be. We obtained Early Science observations 
from a new 4-meter class telescope (Southern Observatory for Astrophysical Research, or 
SOAR) and its imager (SOAR Optical Imaging Camera). We have combined these ob- 
servations with archival observations of the central galaxy in order to compare the spatial 
relationship of the X-ray gas, the radio plasma, the emission-line gas, and the central galaxy 
and to estimate the star formation rate from emission-line and UV indicators. We present 
the details of the observations, data reduction, and results in § |3| including optical imag- 
ing from SOAR and the Hubble Space Telescope 
Double Spectrograph on the Palomar Telescope I 
Array 



3.1), optical long-slit spectra from the 
3.2), radio images from the Very Large 
3.3), X-ray archival data from Chandra (§ 3.4), and ultraviolet archival data from 
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the XMM Telescope (§ 3.5). In § |4l we discuss the implications of these multi-wavelength 



observations, beginning with the relationship of the brightest cluster galaxy and its com- 
panion in § We interpret the morphology of the Ha filaments, partial rims that appear 
to encircle the resolved radio structures in § 4.2, and we discuss further multi-wavelength 



identifications, including the correspondence between the X-ray, optical, and radio in 
the dynamics of the bar and the companion in 



4.3 



4.4[ and the star formation and cooling 

We 



4.5 



rates, including UV measurements from the XMM Optical Monitor camera in 
present a summary of our results in § |5| 

We assume a Hubble constant of 70 km s"^ Mpc^^. At the redshift of 2A0335-I-096, the 
corresponding angular size is 0.7 kpc arcsec"^, and the luminosity distance is 154 Mpc. 



2. The X-ray cluster 2A0335+096 



The X-ray cluster 2A0335+096 {z = 0.035) is a nearby X-ray luminous cluster with a 
cool core and a central radiative cooling time shorter than a Hubble time. The cluster was 



first identified optically (Zwicky et al. 1965). Schwartz et al. (1980) confirmed its positional 



coincidence with an Ariel 1/ X-ray source (Cooke et al. 1978) and noted the short central 



cooling time of the X-ray plasma. Its central galaxy hosts a spectacular, filamentary emission 



line nebula (Romanishin & Hintzen 1988, ; RH88 hereafter). These optical filaments extend 
from a bar-shaped feature in the center of the brightest cluster galaxy (BCG). The BCG 
appears to have a companion galaxy about 4 kpc projected to the NW. Using the ROSAT 



HRI detectors, Sarazin et al. ( 1992 ) revealed a similar bar in the X-rays, and a larger scale 
X-ray elongation of the cluster emission in the same general direction. The BCG is a radio 
galaxy of ~ 33.5 mJy at 1.5 GHz diffused over a region of 29" x 



The BCG nucleus and its companion each have a radio point source (Sarazin et al. 1995). 



(Sarazin et al 


1995 


Sarazin et al. 


1995 


)• 



2A0335-f096 was observed with the Chandra X-ray Observatory (Kawano et al. 2003 



3.4 The Chandra 



Mazzotta et al. 2003). We will discuss these Chandra data further in 
image shows a bar feature, with X-ray deficits to the NE and the SW, identified as possible 
bubbles ( Mazzotta et al.||2003 ). But the details of the X-ray structures and the Ha filaments 
do not coincide, except for the bar around the nucleus, and the extension between the BCG 



and the companion. Mazzotta et al. (2003, MEM hereafter) identified a discontinuity or front 



feature in the X-ray gas to the south of the BCG, confirmed as a cold front by Werner et al. 



(2006). MEM estimated a pressure jump across the feature, and derived a sub-sonic Mach 
number (0.75±0.2), or a propagation speed of about 400(T/2 keV) km s~^. The two peaks in 
the X-ray surface brightness map lie roughly along a line defined by the centers of the BCG 
and its companion. They discussed possible processes explaining the optical filaments in 
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the context of their X-ray observations, ruhng out Rayleigh- Taylor instabihties, but leaving 
open the possibility that the instabilities result from a Kelvin-Helmholtz instability. Such 
instabilities could cause the Ha filaments in a turbulent process. 

Turbulence would result in broadened line profiles and, to some extent, the lack of sys- 
tematic velocity structures across the system. Since the radial velocities of the companion 
and the filaments were unknown to MEM, they could not discuss the dynamics of the inter- 



action. Here, we present in § 3.2 previously unpublished velocity information along the axis 



of the interaction, redshifts and line broadening. 



3. Observations and Data Analysis 

3.1. Narrow-Band iia+ [N II] and Broad-band Imaging 

We obtained broad I-band and narrow-band images (AA ~ 75A, NOAO filter ID 6781- 
78) centered near the redshifted Ha line on Dec. 8, 2005 (UT) with the SOAR Optical Imager 
(SOI), mounted on the SOAR telescope on Cerro Pachon, in Chile. We supplemented these 
data with an R-band (filter F606W) Hubble Space Telescope WFPC2 image. We used the 
WFPC2 associatioiQlD u5a40701b, which is a basic co-addition of two 300-second snapshots 
(Proposal ID 8301, PI Alastair Edge). 

The SOAR data were obtained as part of an Early Science run for the SOI. The pixels 
were binned 2x2, for a pixel scale of 0.154" per pixel. The sky quality was excellent, with 
seeing below 1" and good transparency. The telescope focus needed to be optimized on a 
regular basis (once every two hours) as the object transited. We collected nine 20-minute 
exposures for a total exposure time of 3 hours through the narrow band filter. Each exposure 
was dithered slightly. The narrow-band filter is only 2 inches x 2 inches, so the usable field 
of view is somewhat smaller than the full 5' x 5' FOV of the SOI. After the dithering, our 
useful field of view resulted in an image of 1561 x 1395 image pixels, or 4' x 3.6'. 

A single 10 minute I-band exposure was obtained for the purpose of estimating the 
stellar continuum contribution to the narrow-band observation. All exposures were bias- 
subtracted and flat-fielded with normalized, median-filtered twilight flats taken through the 
appropriate filters. The I-band image exhibited fringing after flat-fielding. We removed 
the fringing signal by subtracting a scaled fringe frame constructed from dark I-band sky 



^This research utilized the faciHties of the Canadian Astronomy Data Centre operated by the National 
Research Council of Canada with the support of the Canadian Space Agency. 
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observations by SOAR observatory personnel. 

We used the IRAF task imalign to align, shift, and trim the images to a common 
image coordinate system using five to six stars in the field as reference points. The shift 
uncertainties were 0.04 pixels. We then combined the narrow-band images by obtaining the 
average in each shifted pixel, after a 3cr clipping filter was applied. This routine removed 
some cosmic-ray features, but we were able to clean single-pixel cosmic-ray events using the 
IRAF task cosmicrays. The pixel scale of the SOI is significantly smaller than the stellar 
PSF, so this process was reasonably successful. 

The world coordinate system (WCS) telemetry from the telescope was not fully im- 
plemented at the time of these observations, so we applied a post-facto astrometric WCS 
solution to the header of the combined narrow-band image using tools available in WC- 
STools (v. 3.6.3, available from the Smithsonian Astrophysical Observatories). We matched 
the positions of 13 objects in the field to a catalog of USN0-A2 stars and galaxies. These 
measurements provided the central astrometric location of the image (J2000), the angular 
scale of 0.154 arcseconds (±0.004"), and a very small angular offset of the field from true 
north. The confidence of the absolute pointing is better than 1", limited by the sparse 
number of actual stars in our small field. Comparisons of features in the archival F606W 
HST/WFPC2 image confirmed the astrometry in the center of the field to better than 0.5". 

Finally, we subtracted a co-aligned, fringe-corrected, scaled (by a factor of 28.57) I-band 
image from the narrow-band image. The scale factor was chosen to subtract the nearby 
cluster galaxies from the narrow band image (Figure [T]). This figure includes indicators 
for the filament structures used to visually define the circles suggested by us in Figure [2} 



discussed further in § 3.3 Since both [N II] and Ho; emission features are included in the 
bandpass of the NOAO filter, we use the long-slit spectroscopy described in §3.2 to apply 
an approximate (~ 15% accuracy) fiux- calibration to the final net emission-line image of 1 
ADU s^^ per 1.26 x 10^^^ erg s~^ cm^^ of pure Ha. This ratio assumes that the mean Ha 
contribution to the total Ha+[N II] emission-line fiux is is 38% throughout the nebula, based 
on the long-slit spectra. 

The emission-line structures are clearly visible in both the full and the the continuum- 
subtracted emission-line images. No filaments or similar structures are visible in the I-band 
image, which is dominated by starlight. To make a direct comparison to fiuxes reported in 
RH88, we measured the total observed Ha emission-line fiux from the nebula within 23" of 
the BCG (masking out a region of r = 20" from the bright star to the southwest) to be 
9 X 10~^^ erg s~^ cm~^. RH88 estimate an emission line fiux of 1 x 10^^^ erg s~^ cm~^, if we 
assume the compatible line ratios. These numbers agree well within our estimated mutual 
absolute calibration accuracies of 10-15%. Enlarging the aperture includes a larger fraction 
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of the detected nebula, so for the purposes of this paper, we will use the a 35" aperture, for 
which we estimate a total pure Ha flux of (1.1 ±0.1) x 10^^'^ erg s~^ cm~^. The uncertainty 
quoted here is based on the uncertain sky background due to sky gradients in the I-band and 
narrow-band frames. An additional uncertainty arises from the scattered light contribution, 
arising because the SOAR telescope was incompletely baffled at the time of the observation. 
Since this scattered-light component is smooth, it does not give rise to spurious structure 
in our image, but it does cause a large-scale gradient across the full field (~ 8% across the 
twilight-fiattened Ha image). We note the total resolved fiux from the nebula is almost a 
factor of 5 larger than that captured in the 2" x 2' spectroscopic slit. 

We note that even though the observation was quite deep, we do not detect filaments at 
even larger projected distances from the BCG. A typical 3a upper limit on such filaments. 



with scales of order 1" — 3" is about 3 x 10 erg s ^ cm 
for total emission line fiux from the Ha-[N II] complex.) 

The Galactic extinction in the direction of 2A0335+096 is Aji 



arcsec ^ of Ha (2.6 times larger 



1.097 (Schlegel et al. 



1998). The total pure Ha luminosity corrected for Galactic extinction but uncorrected 
for any internal extinction, is therefore quite large, ~ 0.8 x 10^^ erg s~^. The total star 
formation rate based on the Ha luminosity, assuming the conversion of SFR(M0 year~^) = 



7.9 X 10~'^^L(Ha) erg s~^ (Kennicutt 1998a) applies in this galaxy, is at least 6 solar masses 
per year. A similar relation for starburst galaxies ( Kennicutt||1998b ) predicts at least 7 solar 
masses per year. 



The lack of calcium emission lines (Donahue & Voit 1993 ) in the nebular gas in 2A0335+096 



indicates that calcium is depleted, and therefore that the emission-line gas itself is dusty. 
This conclusion is nearly independent of the details of the production of the nebular lines. 



Typical dust extinction in these dusty nebula (e.g. Abell 2597; Voit & Donahue (1997)) is 
Ay ~ 1, so correcting for internal extinction would yield a higher star formation rate, up to 
15-20 solar masses a year. 

The B.a + [NII] nebula in 2A0335+096 is similar in luminosity and in appearance to the 
filamentary nebula seen in Perseus's BCG (NGC1275, also 3C84) ( [Fabian et al.][2003| . Some 



of the emission-line structures in NGC1275 have been identified as rising bubbles (Fabian 



et al. 


2003 


Hatch et al. 


2006 



around the radiios source on both sides of the bar feature (Figure |2]). We will discuss these 
emission-line structures further in § 4.2 and the radio contours in § 3.3 The SOAR emission- 
line image is deeper and has better spatial resolution than that of RH88, revealing additional 
emission-line structures south and east of the bar. However, we have detected no emission- 
line counterparts to the X-ray deficits about 38 kpc to the east and 21 kpc to the NW of the 
X-ray peak (MEM). 
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The HST R-band image has several notable features (Figure |3]). The brightest galaxy 
and the companion have smooth elliptical light distributions. The large-scale light distribu- 
tion of the BCG swallows the companion. The stellar light exhibits a distinct, small dust 
lane (1-2 kpc long and < 150 pc wide, barely resolved) north of the companion (Figure [s]). 
This dust lane may extend somewhat closer to the nucleus of the companion than is visible 
in the figure, and a less distinct N-S dust lane, 1" east of the companion nucleus seems 
to be present. These dust lanes appear to be part of the general trail of material behind 
the companion, which suggests that the companion may have had a tiny amount of dust 
in it. The only potential dust feature visible in the BCG is wedge-shaped, south of the 
main nucleus, and pointing in the direction of the nucleus, as indicated in Figure |3| Both 
of these features are located in regions of Ha emission. The wedge feature is filled in with 
high surface-brightness Ha, and the northern dust lane is in the same region as a broader, 
similarly aligned. Ha filament north of the companion. 



3.2. Long-Slit Spectroscopy 

The central galaxy in 2A0335-I-096 was observed by Donahue and Voit using the 5-meter 



Palomar Observatory with the Double Spectrograph (Oke & Gunn 1982) on Dec 31, 1992 
(UT). A composite spectrum of this source was presented by Donahue & Voit (1993). A 2" 
wide slit, approximately 2' in length was placed at a position angle of 140 degrees east of 
north, and centered over the galaxy and its companion. Six 30-minute exposures were taken, 
with the 300-line/mm blue grating centered on 550 nm (2.17 A per pixel) and the 1200- 
line/mm (0.814 A per pixel) red grating centered on 820 nm. The spectral resolution was 
approximately 2A in the red spectrum and 6A in the blue. The angular scale was 0.47" per 
pixel for 24 fim pixels on the red CCD, and 0.78" per pixel for 2 x 15 fim pixels on the blue 
CCD. A composite spectrum from this observation was reported in Donahue & Voit (1993). 
The long-slit image was bias-subtracted and dome-flattened. The wavelength calibration 
was obtained by fitting the known line positions of a neon-argon arc observation and a 
HC observation taken between each sky observation. Star observations taken at different 
positions along the slit and the wavelength calibration were then used to geometrically rectify 
each image. The sky background was removed by identifying night sky regions on either side 
of the object spectrum, then fitting and subtracting a low-order polynomial function from 
the 2D image line by line. A fiux calibration for this instrumental setup was obtained from 
observations of the spectroscopic fiux standards Feige 34 and HD19445, obtained with 6" wide 
slits set to the parallactic angle. The observations were obtained during a 50% illuminated 
moon, but the skies were clear. The seeing was variable and sometimes mediocre (~ 1 — 2"). 



- 9 - 



We have identified four major features along the spectroscopic sht, identified jointly in 
Figure |4j The main nucleus is A, the companion nucleus is B. The filament (S) to the SE and 
the filament (N) to the NW are also labeled. Figure [s] plots the results of fitting each row 
of the 2-D spectrum between the wavelengths of 6750A and 6850A, which includes the Ha- 
[N II] emission line complex. We modeled this complex with three redshifted, identical-width 
Gaussians centered at rest wavelengths of 6548A, 6562. SA, and 6584A, respectively, and a 
flat continuum. The relative emission-line strengths of the two [N II] lines were constrained 
to be in a ratio of 1:3. We fit the spectrum of each image column across the slit. Since the 
pixel scale is 0.47" pixel"^ and the seeing conditions varied from 1 — 2", these plots have 
effectively been passed through a 3-pixel smoothing kernel. The approximate, relative slit 
positions of the image features, as plotted in Figures Figures [5] are reported in the fourth 
column of Table [l] (For reference, arcsecond numbered for Figure [s] corresponds to Line 
63 in the original CCD frame.) 

These features are distinctly seen along the slit as peaks in the Ha line intensity (Fig- 
ure [s]). The locations of these peaks correspond approximately to changes in the velocity 
widths, relative velocity, and [N II] /Ha ratio. The emission line region A', 1.5 arcseconds 
SE of A and surrounding knot B exhibit broader emission lines (FWHM ~ 500 — 600 km 
s~^) than the rest of the nebula (FWHM ~ 200 - 400 km s"^ (Figure [s]) Since both knots 
are associated with radio emission, this broadening may be indicative of optical AGN-type 
activity in both nuclei (both of these nuclei are radio sources). The presence of bright and 
broadened Ha emission is correlated with a noticeable enhancement of the [N II] 6584 / Ha 
ratio in B, but not so much in A and A', where the associated peak in the [N II] ratio 
occurs farther south. (Figure [H]). In the regions where the FWHM > 400 km s~^, the mean 
[N II] /Ha ratio is 1.3 ±0.1 and in the regions where the FWHM was lower, the corresponding 
ratio is 1.1 ± 0.15. An enhanced [N II] line relative to Ha is suggestive of a heating source 



with harder photons than typical of star formation regions (e.g., Baldwin et al. 1981), or 



a possibly a source of supplemental heating above that from photoionization (e.g., Voit & 



Donahue 1997). 



An examination of the relative velocities in Figure |5] shows the presence of at least two 
distinct velocity systems {z = 0.0347), with a velocity difference of about 270 km s~^. The 
recession velocity of Knot A lies (coincidentally perhaps) near the mean, in the center of a 
velocity gradient that extends from A' (at —120 km/s) to B (at +180 km/s). Both A' and N 
are blueshifted about —120 to —150 km s^^ relative to the mean, with the most extreme (and 
faint) velocity component of the NW filament at —180 km s^^. N has somewhat broadened 
emission lines (400 km s~^) with elevated [N II] /Ha as well. Both S and B are receding at 
about +170 km s~^ with respect to the mean velocity of the system. 



The peak of the Ha emission along the sht is somewhat more extended (by a factor of 
two) than a point source. The south end of the peak has the broader emission hues, while 
the north end of the peak lies in the middle of the velocity gradient across the bar. 



3.3. Radio 



We re-reduced archival high-resolution 20-cm A-array data from the Very Large Array 
(VLA), and added it to existing C- and D-array data. This process allowed us to produce 



maps at higher resolution than the 1.5 GHz image presented in Sarazin et al. (1995). We can 
see the inner lobes at a position angle of 65 degrees east of north, oriented perpendicularly 
to the optical emission-line axis angle of 140-145 degrees. An extension of radio emission 
towards the companion is visible along the emission-line bar. The higher resolution radio 
map (Figure [2]) shows the radio-lobes surrounded by curved Ha filaments. 

We estimate the age of the radiating electrons that produce the radio emission following 



Myers & Spangler (1985). We assume that the spectral index between 1.4 and 5 GHz has 
steepened from an initial value of -0.5 to the observed value of about —1 due to synchrotron 
losses (See Table 4 of Sarazin et al. ( 1995[ )). Our results are not very sensitive to whether 
we assume the existence of pitch angle scattering of the relativistic electrons. Assuming the 
magnetic field is at the equipartition value (taking an average value for the two radio lobes 



of 7 jjG, from Table 4 of Sarazin et al. (1995)) we estimate an age of about 25 Myr. On the 
other hand, in some radio sources, there is evidence that the magnetic field is up to a factor 



Carilli et al. 


1994 


Croston et al. 


2005 


Wellman 



et al. 1997). A factor of four weaker magnetic field implies an electron age of about 50 Myr. 



3.4. Chandra X-ray 



We re-analyzed the Chandra data presented by MEM, in order to take advantage of 



the most recent calibrations and to utilize a promising new technique developed by Diehl 



& Statler (2006) for creating X-ray surface brightness and temperature maps. We present 



adaptively binned X-ray and X-ray temperature maps of 2A 0335+096 created from the 



20,000 second Chandra observation (Kawano et al. 2003) (OBSID 919). This observation 



was made with the back-illuminated Advanced CCD Imaging Spectrometer (ACIS) S3 CCD 
on September 6th, 2000. The data were prepared using CIAO v3.3 and the calibration files 



in CALDB v3.2.1. Spectral fitting was performed with XSPEC 11.3.2 (Arnaud 1996) over 
the energy range 0.7 — 7.0 keV. 
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To create the surface brightness and temperature maps, we employed the adaptive 



binning code of Diehl & Statler (2006). Their technique uses weighted Voronoi tessellation 
(WVT) to create spatial bins which maintain a given signal-to-noise (S/N) and uses an input 
background spectrum specific to the observation. The events associated with hot pixels, bad 
pixels, and point sources were not included in the S/N estimation for the definition of the 
spatial bins. Regions of the detector near the edges with less than 40% of the total exposure 
were not included in the final X-ray image. The X-ray surface brightness map was binned 
to a minimum S/N = 5 (Figure [6]). 

The temperature map was extracted from bins defined by an X-ray surface brightness 
map binned to a minimum S/N = 30. The high S/N for the temperature map insures suffi- 
cient counts in each spectral bin to obtain a interesting constraint on the plasma temperature. 
To generate the temperature map, a spectrum was extracted from each WVT-defined bin. A 
corresponding background spectrum for each spatial bin was extracted from the appropriate 
deep background file (provided in the Chandra CALDB by Maxim Markevitcl]|^, which was 
reprojected and matched in gain to the current observation. Each spectrum was fit with a 
projected single-temperature MekaL thermal spectrum and Galactic absorption model us- 
ing a fixed Galactic value of Nh = 1.81 x 10^^ cm~^. The metallicity, temperature, and 
normalization were free parameters. The best-fit temperatures then populate a temperature 
map. The best fit temperatures were also projected onto a uniform grid and interpolated 
using the IDL procedures triangulate and griddata. We confirmed that this interpolation 
did not introduce any spurious features into the temperature map by comparing the same 
contour lines of both maps. The contours of the raw map oscillated around the smoother 
contours of the interpolated map, so we have chosen to show the contours of the interpolated 
temperature map in Figure [7], overlaid on the net Ha image from SOAR. 



3.5. XMM/Optical Monitor: Ultraviolet 



We reduced the X-ray MultiMirror Observatory (XMM) Optical/UV Monitor (OM) 
data of 2A0335+096, taken on August 4, 2003 (DatalD 0147800201, PI Jelle Kaastra). OM 
is a 30 cm telescope mounted on XMM that is capable of simultaneous optical and UV 



observation with the X-ray observation (Mason et al. 2001). Exposures with the UVWl 



(2400-3600 A) and UVW2 (1800-2400 A) filters have been taken. There were five pointings 



http://cxc.harvard.edu/contrib/niaxini/bg/ 
^Interactive Data Language (IDL) is data visualization and image analysis software available from ITT 



Visual Information Systems, http://www.ittvis.coin/ 
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va the UVWl band and eight pointings in the UVW2 band, each with a central small window 
(called "Window 0" with 0.5" pixels) and another bigger one ("Window 1" with 1" pixels 



(Mason et al. 2001). Since the cD galaxy is near the center of the field, it is always present 



in Window of each pointing. The total UVWl exposure is 20600 sec for Window and 
4120 sec for Window 1. The total UVW2 exposure is 34660 sec for Window and 8580 sec 
for Window 1. The cD galaxy 2MASX J03384056+0958119 is a significant UVWl source, 
as well as the companion. There are no UVW2 detections for either source, which is not 
surprising since the telescope sensitivity in the UVW2 band is low and the known Galactic 
extinction in that band and in that direction on the sky is severe. Within the central 7" 
radius (enclosing the light of the companion), the total UVWl luminosity is 5.5 ± 0.4 x 10^^ 
erg s~^, after a correction for coincidence and deadtime loss. We chose to do photometry 
within 7" radius as that aperture encloses most of the UV emission. The level of diffuse UV 
emission beyond the aperture is highly uncertain because of the strong scattered light near 
the middle of the field. The Galactic extinction in the UVWl band has also been applied 



(2.48 mag) based on Cardelli et al. (1989), while the internal extinction is unknown, and 
is set to zero. The contribution from the passive stellar population is subtracted from the 
empirical L^/y^yi - Lj relation derived by Hicks & Mushotzky (2005), resulting in a net 



NUV luminosity excess of 2.9 ± 0.5 x 10 erg s . This excess corresponds to a SFR of 
3.1 — 5.7M0/yr for the assumed IMF with a power-law index of 2.35 — 3.3. 

We have also estimated the UVW2 - UVWl color of > 0.38 mag, without any correction 
for internal extinction. We adopt a Galactic extinction of 3.88 mag at 2050 A, although there 
is quite a range of extinction across the bandpass of UVW2 (3.26—4.14 mag also from Cardelli 
et al. 1989) since the UVW2 band encloses the peak of Galactic extinction curve (~ 2150 



A). This color can be compared with that of Al795's cD (0.0 — 0.15 mag from Mittaz et al. 



(2001). Only a small amount of internal extinction would imply that the intrinsic color of 



2A0335 is similar to that of Al795's cD. 



4. Discussion 



4.1. Merger of the BCG and Companion Galaxy 



Most importantly, from the viewpoint of previous studies of this system, we confirm 
that the companion galaxy to the BCG has a very similar recession velocity (Aw ~ 100 km 



s ^) to the nucleus of the BCG. As noted by Werner et al. (2006) and by Mazzotta et al. 



(2003), the cool X-ray peak, offset from the BCG itself, lies in a line with the BCG and the 



companion galaxy, suggesting that a sub cluster has fallen in along a filament that has fed 
the formation of the cluster. Our radial velocities of the companion and the BCG suggest 



- 13 - 



that the companion and the BCG are not only close in projection but may be in the process 
of merging.. 



4.2. Emission-Line Filament Shapes and Motion of the AGN 

Circles in Figure |2] show the positions and sizes of two apparently circular structures 
visible in the SOAR emission-line image, partially outlined by curved Ha filaments. When 
we complete the circles outlined by these filaments, we obtained the circles on the sky in 
Figure |2] We were struck by the symmetry and the similarity in size of the arcs outlining 
curved structures in the image. The center of the NE circle would be at 3^^ 38™ 41. 15^ 
+09°58'19.9" (J2000). The center of the SW circle would be at 3^^ 38™ 39.8^ +09°58'10.5" 
(J2000). The mid-point of these two positions is 3^^ 38™ 40.5^ +09°58'15.1" (J2000). If the 
AGN is at the peak of the Ha emission, of 3'^ 38™ 12^ +09°58'12" (J2000), then the opening 
angle of the inferred structure, including the nucleus at the origin, is 141.5° in the plane of 
the sky, opening to the northwest, along the bar. 



If the radio plasma in the lobes is 25-50 million years old, as we estimated in § |3.3 the 



separation between the midpoint of the completed Ha arcs and the actual radio/Ha nucleus 
suggests that the AGN is moving through the neighboring X-ray gas at a velocity in the 
plane of the sky of about 55-110 km s~^. This velocity is less than and therefore consistent 
with the 400 km s~^ 3D velocity needed to create the X-ray front feature analyzed by MEM. 
The location of the extended, small-scale, radio emission along what might be the leading 
edge of a cavity supports the suggestion of relative motion between the AGN and the X-ray 
gas. The opening angle between the two lobes of the radio source itself is nearly straight, 
~ 180 degrees. 

While it is possible to draw other circular structures using the filaments as guides, these 
particular structures are the ones closest to and on opposite sides of the AGN. Therefore 
the velocity in the plane of the sky that we infer here is a lower limit. The complexity of 
the emission line nebula appearance indicates that there may have been multiple outburst 
events, occurring along the interaction axis defined by the X-ray cool core, the BCG, and 
the companion. Multiple bubbles and a complex filamentary structure, reflecting a history 
of outbursts from the AGN, have been seen in the X-ray image of M87 in the center of the 



Virgo cluster (Forman et al. 2005) 
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4.3. Other Multi-wavelength Identifications 

The locations of the brightest points in the BCG in the SOAR Ha image, in the 
HST/WFPC2 image, and in the radio image are within 0.5" of each other. We identify 
the nucleus of the BCG with Knot A in the long-slit spectrum (Figure 111 Figure Isl). The 



center of the BCG is also detected in the UV by the Optical Monitor on board XMM (§ 4.5). 
The nucleus of the BCG, as seen by HST, seems to be near the center of the brightest bar 
feature in the emission-line image. The brightest peak in the X-ray emission is south of the 
BCG nucleus by 13.7" or 9.6 kpc. The BCG nucleus corresponds to a fainter peak in the 
X-ray emission. Knot B (Figure |4| |5]) also has counterparts in all datasets, including the 
UV Optical Monitor (XMM), the 3" resolution VLA image, and the Chandra X-ray image. 

The X-ray emission does not show a one-to-one structural correspondence with the Ha 
emission (Figure [6]), but there is a choppy ridge of high X-ray surface brightness that follows 
the more pronounced Ha bridge between the two galaxies. Just outside the Ha nebula, to 
the south, and past the X-ray peak, there is an X-ray cool front (MEM). It's impossible to 
see the behavior of the optical emission line gas in this region in the SOAR image because of 
a very bright star, but even the Palomar spectrum, which is less affected by the bright star, 
does not show much hint of extended emission in that direction. Beyond the Ha nebula, 
the X-ray surface brightness profile drops off suddenly, as if a minimum of X-ray surface 
brightness were needed to host Ha filaments. 

Our temperature map of the cluster shows that the Ha emission-line region lies entirely 
within the 2.5 keV contour (Figure [?]). This X-ray gas has the lowest entropy and shortest 
cooling time in the cluster, and it has the same elongated morphology as the Ha nebula. 
This correspondence suggests that the cool X-ray core and the BCG are related, even though 
the X-ray peak is not at the center of the BCG. We note that the X-ray surface brightness 
peaks in two cool (~ 1.5 keV) clumps to the SW, one which may be associated with the SW 
filament and the other where the foreground star obscures our view of the emission-line gas. 
Thus the coolest X-ray gas does not correspond exactly to the Ha nebula, a fact we will 



return to in § 4.5 As mentioned in § 3.3 and in § |4.2[ the radio source may be affecting the 



optical filaments on either side of the optical/radio nucleus, as seen in the Ha image. From 



this analysis, we concluded in § 4.2 that the AGN may be in motion with respect to the 
optical nebula. This hypothesis is supported by analysis of the emission-line spectra, which 
we discuss next. 
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4.4. Dynamical Clues: Systematic Velocities along the Bar and Between the 

Knots 



We discuss here possible interpretations of the systematic velocity trends we have mea- 
sured along the bar of optical line emission and surrounding filaments: a rotating filament 
system in a counter-rotating disk or an interacting system, with a bar of stripped gas. We 
also discuss the relationship to the motion of the cool clump seen in the X-ray images. 



This bar has been interpreted as an edge-on disk by both Romanishin &; Hintzen ( 1988 ) 



and Sarazin et al. (1995). Sarazin et al. (1995) suggest that this disk may determine the 



orientation of the radio source. Romanishin & Hintzen (1988) suggest a compact 2 kpc 



disk inside a 17 kpc bar. Neither of these groups had spectroscopic information about the 
velocities along this structure, however. It is possible to interpret the velocity structure 
along the bar as a large, rotating filament system around Knot A with a circular velocity 
of 170 km s~^ (using the extreme ends of the bar N and S to define the "rotation"), and a 
diameter of about 20" or 14 kpc. This circular velocity and size implies a gravitating mass 
of at least 5 x lO^'' M0. The center of this system lies close to the center of the BCG. 



The inner structure of the bar could be a counter-rotating disk as suggested by Roman 



ishin & Hintzen (1988), with a diameter of about 7.5" or 5.25 kpc and a circular velocity 



of ~ 120 — 150 km s~^. The mean redshift of the bar may be somewhat greater than that 
of the more extended system, but only by about 30 km s~^. The intensity peak of the 
Ha emission (Knot A) lies at or near the center of this bar. The implied mass interior to 
r = 3.5 kpc is approximately 1.4 x 10^" Mq. A similar, but larger-scale structure was seen 



in 3C275.1 (Hintzen & Stocke 1986), but we note that even the gradient in this structure 



was not considered strong indication of rotation by these authors. 

A second, possibly more compelling interpretation of the velocity system is that it is the 
result of an stripping interaction between the two nuclei. Knots A and B. Stripping tails have 



been seen in the X-ray by Sun et al. (2006) in Abell 3627 and in the optical (Gavazzi et al. 



2001) in Abell 1367; see also Sun et al. (in preparation) in Abell 3627. The portion of the 



longslit spectrum with the broadest emission lines seems to lie closer to the southern end of 
this bar, not in the center of it. The bar of emission between Knots A and B may be stripped 
material, with a range of velocities along the stream, stretching from Knot A towards B. If 
the overall direction of motion is from NW to SE, as suggested by the orientation of the cool 
front and the BCG, it is possible that the BCG, moving in the same general direction as the 
cool core, passed through B, stripping material from B in the process. 

A tidal origin for the bar seems unlikely due to the lack of a similar feature in the 
stellar continuum image, unless the tides also induce physical processes that light up the 
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emission-line and X-ray gas without significantly affecting the I-band light distribution from 
stars. 

Knot A and the companion galaxy (Knot B) have a velocity separation of about 100 km 
s~^ along the line of sight, and a projected separation of 6.6" or 4.6 kpc. The crossing time 
therefore is several times 10'' years. (We have made the simplifying approximation that the 
physical separation of the two knots is -\/3 larger than the projected separation, and that 
there is no motion in the plane of the sky.) The gravitational free fall time {tjf ~ (Gp)^^^^) 
for a system with M ~ 10^° Mq and r = 7 kpc is about 2 x 10^ years. 

Interesting dynamical information also comes from the X-ray image. Recall the shape 
of the X-ray edge and pressure differential imply that a cool clump is moving to the south 
with a Mach number of ~ 0.75 (MEM). This Mach number implies that the cool clump is 
moving with respect to the larger cluster at ~ 500 km s~^. The relative position of the 
circular Ha filaments and the AGN suggest that the AGN may be moving with respect to 
the local ICM in the same direction as the cool clump, in the plane of the sky. The filaments 
defining the circular structure are arranged symmetrically around the radio source, with a 
slight offset suggesting lower limit to the velocity in the plane of the sky of ~ 50 — 100 km 
s~^. These clues indicate that the AGN may be moving in the same general direction as the 
cool clump, but possibly at a lower speed. 



4.5. Star Formation Rates and Mass Cooling Rates 

One problem with the original cooling flow model was the large discrepancy between 



the high X-ray mass cooling rates and the low star formation rates (SFRs). In § 3.1 



we 



derived a minimum SFR (with no intrinsic dust extinction correction) of at least 7 M0 yr ^ 



from the SOAR Ha imaging data. Peterson et al. (2003) constrained the cooling rate of the 



X-ray gas in the core of 2A0335+096 with measurements and limits on X-ray emission lines 
seen in the spectrosopic grating data from the XMM RGS instrument. Individual X-ray 
lines provide independent constraints on the cooling rate, so the cooling constraints depend 
on the spectral regime under consideration. They derive an M(0.8 — 0.4 keV) = 20±10 Mq 
yr~^ and M(0.4 — 0.2 keV) < 84 Mq yr^^. While the uncertainty is large, this updated 
cooling rate of the X-ray gas is close to the SFR in the cD galaxy. 

From the UV observations, we derived a near UV- luminosity excess of 2.9 ± 0.5 x lO''^ 
erg s~^. This excess corresponds to a SFR of 3.1 — 5.7M0/yr for the assumed IMF with a 
power-law index of 2.35 — 3.3 and zero internal extinction. If we apply an internal extinction 



correction corresponding to Ay ~ 1.0 (typical for cooling flow nebulae, (e.g., Voit & Donahue 
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1997)), the estimated SFR would be up to 10 times larger. The rate inferred without internal 
extinction is somewhat less than that inferred from the Ha emission, but that is to be 
expected since even tiny amounts of extinction would increase the UV estimate substantially. 
We note the Ha emission spans a much larger region (r ~ 35" in radius) than that of the UV 
emission (within 7" radius), but that most of the Ha is concentrated close in. We conclude 
that the observed UV excess is consistent with the SFR inferred from the Ha observations. 

Yet another independent constraint on star formation was sought from the lUE data 
archive. We inspected two archival lUE observations of 2A0335+096, SWP24934 and SWP43531. 
The short wavelength (SWP) camera on lUE, a 45-cm UV telescope (1978-1996) was sen- 
sitive to UV emission from 115-200.0 nm in a large aperture of 10" x 20", blueward of the 
XMM/OM UVWl band. An upper limit to UV continuum was reported from the former 



observation of 10,800 seconds by Crawford & Fabian (1993). We analyzed a longer, later 



lUE observation, of 24,000 seconds, but we make only a marginal improvement on Craw- 
ford et al.'s original upper limit. Cosmic-ray events obscure the spectrum at the position of 
redshifted Lyman-a at the nominal pointing position. A 3a limit of < 2.4 x 10~^^ erg s^^ 
cm~^ A~^, measured between 128 and 136 nm, was obtained, corresponding to a Galactic 
extinction-corrected flux limit of 6 x 10"^"^ erg s~^ cm~^ A~^. (No emission, extended or 
otherwise, is visible in the lUE field.) A Lyman a line with a width of 500 km s~^ would be 
unresolved in an lUE spectrum (1.67 Angstroms per element), so a 3a line would have to be 
about 4 X 10~^^ erg s~^ cm~^, corresponding to an upper limit of about 5 x 10^^ erg s~^ on 
Lyman a luminosity and < 1.3 x 10^^ erg s~^ on continuum emission in a bandpass between 



128-136 nm, corrected for an extinction of E{B — V)= 0.41 and a Cardelli et al. (1989) law. 
This limit is about 50% of what one would expect based on the Ha luminosity we estimate 
for the system with no internal extinction. However, the amount of dust required to push 
the Lyman-a line below detectability by these lUE data is extremely small. Therefore, the 
lack of detection of either a Lyman-a line or an ultraviolet continuum from 2A0335-I-096 by 
lUE, together with detections by the XMM OM and SOAR, indicates that there is some 
internal dust in the BCG of 2A0335+096. 



The total amount of molecular gas in the galaxy was estimated by Edge & Frayer (2003 ) 
based on CO observations to be 1.1 ± 0.4 x IO^/i^q Mq. Such a gas supply, if not replenished 
by cooling or accretion, would be completely used up by star formation in 100 million years. 
Therefore, this system could be classified as a starburst galaxy. However, since the star 
formation rate is similar to the estimated cooling rate from the hot ICM by Petersen et al. 
(2003), we suggest that it is possible the gas supply is being replenished by cooling. 

Radio emission is thought to be related to accretion onto an AGN. This accretion may 
be fueled by a recent interaction with a gas-rich system. It is possible that a galaxy-galaxy 
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interaction has induced star formation in gas accumulated from the ICM or an ISM gas 
source. We point out that a scenario where hot ICM has cooled and started to form stars is 
less tenable if the star-forming gas turns out to be dusty (Donahue & Voit 1993). Since the 
dust sputtering time in hot gas from Draine & Salpeter (1979), the hot gas in the core of a 
typical cluster sputters refractory grains with size a < 0.1/im (the majority of dust grains 
by number) in about 2 x lO'^na ^ao.i years where n2 = njy/10~^ cm~^ and Oo.i = a/O.l/xm, 
so theoretical prejudice is that such cooled ICM gas would be dust-free. 



5. Summary and Conclusions 

We summarize our results as follows. 



1. Our SOAR Ha image of 2A0335+096 is deeper, more sensitive, and has better seeing 
than the one presented in RH88. We clearly detect the emission-line filaments on the 
both the north and south sides of the galaxy, and we see filamentary details in excellent 
relief. We don't detect much more Ha beyond ~ 25 kpc from the central galaxy. These 
observations limit the presence of somewhat lower low-surface brightness emission-line 
gas at large radii (~ 3 x 10^^'' erg s~^ cm~^ arcsec"^ on scales of 1" - 3".) 

2. We present velocity information along the interaction axis defined by the X-ray cool 
front, the BCG, the X-ray and Ha bar, and the companion galaxy that show the 
companion galaxy is likely to have interacted with the BCG. The velocity structure 
along the X-ray and optical emission-line bar suggests that ram pressure stripping 
might be a relevant process here. 

3. In the frame of the AGN (indicated by the center of the BCG, Knot A), the N/NW 
filament /blob is moving toward us, and the S/SE filament /blob is moving away from 
us. However, a 2-D velocity map would be required to confirm whether the large-scale 
nebula is rotating in the frame of the AGN. We cannot tell whether the large-scale 
system is rotating with these data. 

4. The main evidence for a earlier interaction between A and B comes from the morphol- 
ogy of the emission bar extending between A and B, and the continuous transition of 
radial velocities from A to B seen in the long-slit spectrum. 

5. Because the companion (Knot B) seems to have already interacted with the BCG, we 
conclude that Knot B is now beyond the BCG, pulling out a bar of emission-line gas 
with it. Since both galaxy cores are also emission-line sources, it is impossible to say 
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whether B is stripping gas from A, vice- versa, or both galaxies are stripping gas out of 
each other. The passage of the companion galaxy through the BCG may have induced 
star formation possibly in pre-existing gas accumulated from a cooling flow or in cold 
gas the companion may have contributed to the system. A stream of emission-line 
gas and dust extends past the companion to the north, suggesting that some of the 
B companion's motion is also in the plane of the sky, and B may be now proceeding 
back towards the SE, towards A. One hint that this may be the case is that the gas 
in the NW filament has the same systematic velocity as A. Possibly ICM interactions 
and gravity has decelerated that gas trailing companion B, and so its relative velocity 
to the BCG is now zero. 

6. The interaction may have also triggered a radio source by providing fuel or by exciting 
instabilities in gas orbiting the central supermassive black hole; however, this scenario 
cannot be distinguished from fueling the radio source with gas condensing from the 
ICM. The higher- resolution radio image we present here shows two close-in lobes of 
radio plasma. The optical cmission-linc filaments filaments near the lobes seem to arc 
around the radio plasma. The higher resolution radio data from the VLA and the 
deeper, higher resolution emission line imaging from SOAR provide complementary 
insight to this source. 

7. We have derived a star-formation rate from the UV/OM images of 2A0335-I-096. These 
rates are marginally consistent with, and lower than, constraints on the X-ray mass 
cooling rate from XMM RGS spectra. The approximate total Ha luminosity is also 
consistent with the estimated star formation rate and rate of gas cooling from X-ray 
emitting temperatures. This consistency suggests in this system at least that the star 
formation is a possible sink for gas cooling from X-ray temperatures. 

We note an interesting similarity in three timescales: (1) the consumption or production 
timescale of the molecular gas (the quantity of the cold gas divided either by the X-ray cooling 
rate or the star formation rate), (2) the estimated age of the radio source derived from the 
synchrotron age of the electrons, and (3) the approximate time of last interaction between the 
companion and the BCG. The similarity in timescales suggests that the common triggering 
mechanism for the the radio source and the starburst in this system may be the interaction 
of a companion with a BCG residing in low-entropy ICM gas. 

Our detailed study of this cluster core reveals details about this system that suggest 
the interaction of a companion galaxy with the BCG is the trigger for the radio source and 
the starburst event. Whether or not the cooling fiow was required to supply the gas that 
feeds the black hole and that supplies the starburst is uncertain, but the mass cooling rates 
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inferred from the XMM RGS spectrum are consistent with the star formation rates inferred 
from Ha and UV observations presented here. A two-dimensional study of the emission hne 
system, possibly with an integral field unit, could further distinguish the dynamical models 
suggested by our Ha data. 
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SOAR Ha+[N II] 




Fig. 1. — A logarithmically-scaled version of the continuum-subtracted emission-line image 
of 2A0335-I-096 from SOAR Optical Imager. North is up and East is to the left in this 
figure. The nuclei of the BCG and the companion are labeled A and B respectively. The 
emission-line image shows both emission-line knots embedded in filaments. The filaments 
curve in arcs that extend above and below the emission line bar. The "arcs" indicated here 
are higher surface brightness filaments surrounding the newly resolved radio source, and are 
used to define the circles suggested in Figure [2] 
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Fig. 2. — Radio contours of the highest resolution VLA image (3") overlaid on the optical 
emission line greyscale image. The radio data were obtained Sept 22, 1987 and reprocessed 
on March 30, 2006. The (red) radio contour levels are 0.5, 1.0, 1.5, 2.0, 6.0 mJy per beam. 
The black circles are the locations of the putative circular emission-line features described in 
the Discussion, 8.25" and 8.29" in radius. The NE structure is on the left, the SW structure 
is on the right. The peak of the radio source aligns within 1" of the peak of the emission-line 
source in the BCG. 
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Fig. 3.— A grey scale image of 2A0335+096 from the archival HST WFPC2 red (F606W) 
filter imaging. North is up and East is to the left in this figure. The grey scale image to 
the right is a high contrast scaled image of the same image, indicated by the box on the left 
image. 
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Fig. 4. — The position angle and slit position position for the next 4 plots. The line on 
the 2-D long-slit spectrum is 60 pixels long, starting from the pixel on the lower left hand 
side. The length of that line matches the length of the slit pictured on the right-hand image 
showing the angle and width of the spectroscopic slit. The actual full slit was 2 arcminutes 
in length and 2 arcseconds wide. The 60 pixels (29 arcseconds) spanning the length of the 
slit was used for the flux calibration of the SOAR emission-line image. 



- 29 - 



SxlO 



6x10 



X 4x10 
r3 



2x10 



-16 



-16 



-16 



-16 





S A 




A 


H 


N 




\ 












/ 











10 20 
Arcseconds 



50 



X 




10 20 
Arcseconds 



50 



600 



400 



200 



u 
_o 

> 




10 20 
Arcseconds 



50 



o 
_o 

> 
'o 

!^ 
-(— ' 
C 

(D 
O 




10 20 
Arcseconds 



50 



Fig. 5. — Upper Right. Ha intensity as a function of slit position. The brightest peak of 
the Ha emission is centered in the BCG stellar continuum. The locations of the knots are 
marked by vertical lines in all the figures. Knots A and B are the centers of the BCG and 
the companion. Knot A' is the location of the gas with the broadest line profiles, located 
somewhat SE of the centroid of the brightest Ha peak. N and S mark the location of the 
NW and SE gas clouds in the sht. Lower Right. The mean width of the Ha and [N 11] 
lines in km s"^ at each position along the slit. A Gaussian shape was assumed to estimate 
the FWHM from the width a of the Gaussian lines, where FWHM = 2(2 log 2) VV. Upper 
Left. The ratio of the integrated [N II] 6584 line to the Ha line, assuming that the lines were 
the same width. Each line of the spectrum along the slit was fit separately. Seeing correlates 
the measurements across 2-3 pixels. Three distinct peaks are seen. Two are associated with 
the brightest peaks, and the third is associated with the NW filament. Lower Left. The 
mean recession velocity of all 62 positions along the slit was subtracted from the best-fit 
velocity to obtain a relative velocity for each slit position, in km s~^. The mean redshift 
measured from these spectra and used here is 0.0347. 
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Fig. 6. — Left, a color-coded image of the adaptively-binned (Sex) X-ray (0.7-7.0 keV) surface 
brightness distribution for 2A0335-I-096. The color bar corresponds to the mean photon 
counts per (0.5" x 0.5") pixel in the X-ray image. The net exposure time was about 20,000 
seconds. Note that the nuclei of the galaxies A and B are not located at the brightest X-ray 
peak, but are associated with fainter X-ray peaks. Right, a contour map of the X-ray surface 
brightness isophotes are plotted over the net Ha image from SOAR, at the same scale and 
orientation (North up. East left) as the left image. While the bar of Ha surrounding the BCG 
nucleus does appear to have a counterpart in the X-ray emission, there is no significant one- 
to-one correlation of X-ray and Ha surface brightness. We note that the Ha appears to be 
limited to the region where the X-ray surface brightness is greater than about 3 photons per 
pixel, which is an approximate X-ray surface brightness of 3 x 10~^^ erg s~^ cm~^ arcsec"^. 
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Fig. 7. — A grey-scale image of the net emission-line image, with the contours of the inter- 
polated temperature map superimposed. The temperature map contour levels are 1.75, 2.0, 
2.5, and 3.0 keV, inner to outer (the contours map a decrement in X-ray temperatures). The 
Ha filaments appear to be confined to a region where the X-ray temperature is less than 
2.5 keV, and the coolest lumps (the innermost contours, where < 1.75 keV) correspond to 
the brightest Ha regions, except for the cool X-ray lump obscured by the bright star to the 
south-east, which cannot be seen optically. North is up. East is left. 
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Table 1. Summary of Emission-Line Features 



Label 


RA (J2000) 


Dec (J2000) 


Mean Slit. Pos. 


Redshift 








(arcsec) 






A 


03:38:40.6 


+09:58:12 


10.0 


0.0347 (0.0343-0.0351) 


Bridge A-B 








0.0352 (0.0352-0.0353) 


B 


03:38:40.3 


+09:58:18 


16 


0.0353 (0.0351-0.0353) 


NW Fil.^ 


03:38:40 


+09:58:22 


22.5 


0.0343 (0.0341- 0.0347) 


SE Fil.b 


03:38:41 


+09:58:04 


2.0 


0.0352 (0.0351-0.0353) 



^Each pixel in the spatial direction is 0.47 arcseconds. The slit width was 2 arcseconds. 
'^These astrometric locations are approximate, within 3". 



